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We present measurements of branching fractions and charge asymmetries in _B-meson decays to 
p"'"7r°, p^tt"*" and p'\^ . The data sample comprises 89 x 10^ T(4S) — » BB decays collected with the 
BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We find the charge-averaged 
branching fractions B{B+ p+7r°) = (10.9 ± 1.9(stat) ± 1.9(syst)) x 10"^ and B{B+ p°-K+) = 
(9.5 ± 1.1 ± 0.8) X 10-^ and we set a 90% confidence- level upper limit Z3(B° p"tv°) < 2.9 x lO"'*. 

We measure the charge asymmetries A^p"" = 0.24 ± 0.16 ± 0.06 and A^°p^ = -0.19 ± 0.11 ± 0.02. 
PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The study of B-meson decays into charmless hadronic 
final states plays an important role in the understanding 
of CP violation in the B system. Recently, the BABAR 
experiment performed a search for CP-violating asym- 
metries in neutral B decays to p^n^ final states 0, 
where the mixing-induced CP asymmetry is related to 
the angle a = arg [-Vt^Vtl/V^^V;,,] of the Unitarity Tri- 
angle The extraction of a from p^n^ is complicated 
by the interference of decay amplitudes with differing 
weak and strong phases. One strategy to overcome this 
problem is to perform an SU(2) analysis that uses all 
pn final states |3j. Assuming isospin symmetry, the an- 
gle a can be determined free of hadronic uncertainties 
from a pentagon relation formed in the complex plane 
by the five decay amplitudes P^t^^, B^ /0^7r+, 

5° ^ p07r°, B+ p+ttO and B+ p°Tr+ ||. These 
amplitudes can be determined from measurements of the 
corresponding decay rates and CP-asymmetries. The 
branching fractions have been measured for —> p~^tt~ 
and B^ p^n^, and an upper limit has been set for 

B"^p"^° Hi. 

In this letter we present measurements of the branching 
fractions of the decay modes B~^ p+7r° and i?+ 
p^TT^, and a search for the decay p°7r°. All three 

analyses follow a quasi- two-body approach 0,0]. For the 
charged modes we also measure the charge asymmetry, 
defined as 

(1) 

T{B- ^ /) + r(B+ ^ /) ^ ' 

where / and / are the final state and its charge- 
conjugate, respectively. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy 
e+e" storage ring at SLAC. The sample consists of 
(88.9 ± 1.0) X lO^ BB pairs collected at the r{AS) reso- 
nance ( "on-resonance" ) , and an integrated luminosity of 
9.6 fb"^ collected about 40 MeV below the TiAS) ("off- 
resonance" ) . 

Each signal B candidate is reconstructed from three- 
pion final states that must be Tr+Tr^Tr", 7r"'"7r~7r"'", or 
7r'''7r~7r''. Charged tracks must have ionization-energy 
loss and Cherenkov-angle signatures inconsistent with 
those expected for electrons, kaons, protons, or muons 0. 
The TT*^ candidate must have a mass that satisfies 0.11 < 
to(77) < 0.16GeV/c^, where each photon is required to 
have an energy greater than 50 MeV in the laboratory 



frame and to exhibit a lateral profile of energy deposi- 
tion in the electromagnetic calorimeter consistent with 
an electromagnetic shower 0- The mass of the recon- 
structed p candidate must satisfy 0.4 < r7i(7r+7r°) < 
1.3GeV/c2 for p+ and 0.53 < m{TT+TT~) < O.gCeV/c^ 
for p". The tight upper m(7r+7r^) cut at 0.9GeV/c^ 
is to remove contributions from the scalar /o(980) res- 
onance, and the tight lower cut is to reduce the con- 
tamination from ifg decays. To reduce contributions 
from B'^ — > p+TT^ decays, a ^ p°7r° candidate 
is rejected ff 0.4 < m(7r±7r°) < 1.3GeV/c^ For the 
£?+ p+tt'' and B'^ p^n^ modes, the invariant mass 
of any charged track in the event and the tt*^ must be 
less than 5.14GeV/c^ to reject B'^ tt+tt'^ background. 
For the B^ — > p'^Tr^ mode, we remove background from 
charmed decays B D'^X, K~^w~ or tt+tt", by 

requiring the masses m(7r+7r~) and m{K^TT~) to be less 
than 1.844 GeV/c^ or greater than 1.884 GeV/c^. We take 
advantage of the helicity structure of B pn decays by 
requiring that | cos^pj > 0.25, where 6p is the angle be- 
tween the tt" (tt^) momentum from the p^ (p") decay 
and the B momentum in the p rest frame. 

Two kinematic variables, AE and mss, allow the dis- 
crimination of signal B decays from random combina- 
tions of tracks and tt" candidates. The energy difference, 
AE', is the difference between the e~^e~ center-of-mass 
(CM) energy of the B candidate and \/s/2, where ^/s 
is the total CM energy. The beam-energy-substituted 
mass, iriES, is defined by y/{s/2 + pi ■ pB)^/Ef - p|, 
where the B momentum, p^, and the four- momentum 
of the initial state (i?j, Pi) are measured in the labora- 
tory frame. For B~^ p°7r+ we require that —0.05 < 
AE < 0.05 GeV while for both modes containing a 7r° 
we relax this requirement to —0.15 < AE < 0.10 GeV. 
For both 5+ p'^7r+ and 5° — > p°7r° we require that 
5.23 < ruEs < 5.29GeV/c2 while for B+ p+n° it is 
relaxed to 5.20 < mijs < 5.29 GeV/c^ 

Continuum e+e~ ~^ QQ il = u,d,s,c) events are the 
dominant background. To enhance discrimination be- 
tween signal and continuum, we use neural networks 
(NN) to combine six discriminating variables: the re- 
constructed p mass, I cos^pl, the cosine of the angle be- 
tween the B momentum and the beam direction in the 
CM frame, the cosine of the angle between the B thrust 
axis and the beam direction in the CM frame, and the 
two event-shape variables that are used in the Fisher dis- 
criminant of Ref. 8] . The event shape variables are sums 
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TABLE I: Numbers of selected events from on-resonance data, 
signal efficiencies, relative fraction of misreconstructed and 
wrong charge events from MC. 





B+ ^ p+7T° 


B+ pV+ 


B° ^ pV" 


Selected events 


13f77 


855f 


7048 


Signal efficiency 


f7.5±0.f% 


28.3 ± 0.1% 


20.0 ±0.1% 


Misreconstructed 


38.6 ± 0.2% 


7.1 ±0.1% 


9.1 ±0.2% 


Wrong charge 


8.f ±0.f% 


1.6 ±0.1% 





over all particles i oipiX 10056*^1", where n = or 2 and 0i 
is the angle between momentum i and the B thrust axis. 
The NN for each analysis weighs the discriminating vari- 
ables differently, according to training on off-resonance 
data and the relevant Monte Carlo (MC) simulated sig- 
nal events. The final pir candidate samples are selected 
with cuts on the corresponding NN outputs. 

To further discriminate further between signal and con- 
tinuum background, for the B'^ — > p°7r" mode, we use the 
separation between the vertex of the reconstructed B and 
the vertex reconstructed for the remaining tracks. This 
separation is related to At, the difference between the 
two decay times, by Az = c/3'^At, where for PEP-II the 
boost is (3j — 0.56. 

Approximately 33%, 7%, and 8% of the events have 
more than one candidate satisfying the selection in the 
B+ p+7r°, B+ p°'!T+, and B° p^ir^ decay mode, 
respectively. In such cases we choose the candidate with 
the reconstructed p mass closest to the nominal value of 
O.TTGeV/c^. Table d summarizes the numbers of events 
selected from the data sample and the signal efficien- 
cies estimated from MC simulation. Some of the actual 
signal events are misreconstructed; this is primarily due 
to the presence of random combinations involving low 
momentum pions. For the charged B modes we distin- 
guish misreconstructed signal events with correct charge 
assignment from those with incorrect charge assignment. 
These numbers, estimated from MC, are also listed in 
Table H 

We use MC-simulated events to study the background 
from other B decays, (B-background), which include 
both charmed (b c) and charmless decays. In the se- 
lected p+7r° {p°Tr+, p°Tr°) sample we expect 205±46 (73± 
19, 59 ±18) fe^ cand 228 ±77 (92 ±11, 74 ±22) charm- 
less background events. All the three analyses share the 
major i?-background modes: B'^ — > p+Tr^, longitudi- 
nally polarized B'~' P~^P~, and i?+ P^P^- Other 
important modes include _B+ p+vr*' (for B'^ P^'^'^), 
B+ (ai7r)+ (for B+ p+tt"), B+ K*{892fTr+ (for 
_B+ —> p'^7r+), and background modes containing higher 
kaon resonances. 

An unbinned maximum likelihood fit is used for each 
analysis to determine event yields and charge asymme- 
tries. To enhance discrimination between signal and 



background events, we use the B-flavor-tagging algo- 
rithm developed for the BABAR measurement of the CP- 
violating amplitude sin 2/3 8], where events are sepa- 
rated into categories based on the topology of the event 
and the probability of misassigning the _B-meson flavor. 
The likelihood for the Nk candidates tagged in category 
k is 

Nk f Nb ^ 

where N^^ is the number of signal events in the en- 
tire sample, is the fraction of signal events tagged 
in category fc, iV^'' is the number of continuum back- 
ground events that are tagged in category fc, and Nb is 
the number of _B-background modes. N'^ is the sum of 
the expected event yields for signal {etN'''"), continuum 
(iV^') and fixed B background. For the charged modes 
the asymmetries are introduced by multiplying the sig- 
nal yields by 5(1 — QiAcp), where Qi is the charge of B- 
candidate i. The likelihood term Cf- ^ corresponds to the 
jth i?-background contribution of the Nb i?-background 
classes. The total likelihood is the product of likelihoods 
for each tagging category. 

The probability density functions (PDF) for signal 
and continuum, V^^ and "Pf, are the products of the 
PDFs of the discriminating variables. The signal PDFs 

are given by v\!''"^^ = V^p''''^ (mEs) ■ V^p'^'^^ {AE) ■ 
■p^p^r) (]\[fs[) for the charged B decay modes, and by 
pp'V ^ Vp"^\mEs)-VP°''\AE)-v(''\'^^)-vf''\At) 
for _B" — > p^tt'^. Each signal PDF is decomposed into two 
parts with distinct distributions: signal events that are 
correctly reconstructed and signal events that are misre- 
constructed. For the charged B modes, each PDF for the 
misreconstructed events is further divided into a right- 
charge and wrong-charge part. The niESi and NN 
PDFs for signal and for B background are taken from 
MC simulation. For continuum, the yields and PDF pa- 
rameters are determined simultaneously in the fit to on- 
resonance data. 

In the B^ p°7r° decay the At distributions for signal 
and B background are modeled from fully reconstructed 
B^ decays from data control samples 8] . The continuum 
At parameters are free in the fit to on-resonance data. 

To validate the fit procedure, we perform fits on large 
MC samples that contain the measured number of signal 
and continuum events and the expected _B-background. 
Biases observed in these tests are largely due to corre- 
lations between the discriminating variables, which are 
not accounted for in the PDFs. For p+Tr*^ and p^ir^ they 
are not negligible and are used to correct the fitted sig- 
nal yields. In addition, the full fit biases are assigned as 
systematic uncertainties on all three signal yields. 

Contributions to the systematic errors are summarized 
in Table ^ Uncertainties in the signal MC simulation 
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TABLE II: Summary of the systematic uncertainties. 



Error source 




(events) 


pV 


^CP 

(10 


^CP 


Signal model 


10.7 


3.8 


3.3 


3.4 


0.3 


Fit procedure bias 


14.4 


8.2 


2.0 






B background 


11.2 


2.3 


3.3 


5.0 


2.2 


Detector charge bias 








1.0 


0.9 


Total fit error 


21.1 


9.3 


5.1 


6.1 


2.4 


Relative efficiency error 


11.6% 


7.2% 


7.0% 







are obtained from a topologically similar control sam- 
ple of fully reconstructed — > D~ decays. For the 
i3+ — > p+7r° channel we also use i?+ K^n^ decays 
to estimate the uncertainty in the Ai? model. We vary 
the signal parameters, that are fixed in the fit, within 
their estimated errors and assign the effects on the sig- 
nal yields and charge asymmetries as systematic errors. 
The expected yields from the _B-background modes are 
varied according to the uncertainties in the measured 
or estimated branching fractions. Since _B-background 
modes may exhibit direct CP violation, the correspond- 
ing charge asymmetries are varied within their physical 
ranges. For — > p^tt*^, the systematic uncertainty due 
to interference with —> p^Tr~ is found to be 1.5 events. 
This is obtained by repeating the fit to data, after remov- 
ing the cut on r7i(7r^7r*'). Systematic errors due to possi- 
ble nonresonant B^ — > 7r+7r"7r*' decays are derived from 
experimental limits 0. Contributions from nonresonant 



_B+ vr+vr'^Tr" for the p+vr" mode and _B+ tt+tt 7r+ 



for the p°7r+ mode are estimated to be negligible. For the 
i3+ /9°7r+ and i?° p^n° decay modes, systematic 
uncertainties due to interference between and /o(980) 
or a possible broad scalar cr(400 — 1200) were also stud- 
ied and found to be negligible. Repeating the selection 
and fit for all three modes, without using the p-candidate 
mass and helicity angle, gives results that are compatible 
with those reported here. In the B^ p°7r"'" case, the 
analysis was repeated in the region |cos0p| < 0.25, and 
the resulting signal yield was consistent with zero. 

After correcting for the fit biases we find from the max- 
imum likehhood fits the event yields, N{p+7r") = 169.0± 
28.7, iV(p''7r+) = 237.9±26.5, and N{p°7r") = 24.9±11.5, 
where the errors are statistical only. Figure ^ shows dis- 
tributions oi niES and enhanced in signal content 
by cuts on the signal-to-continuum likelihood ratios of 
the other discriminating variables. The statistical signif- 
icance of the previously unobserved B^ p+Tr" signal 
amounts to 7.3ct, computed as v^SSIogZ, where Alog£ 
is the log-likelihood difference between a signal hypoth- 
esis corresponding to the bias-corrected yield and a sig- 
nal hypothesis corresponding to a yield that equals one 
standard deviation of the systematic error. We find the 




to 




FIG. 



in p^TT 



1: Distributions of niES and AE for samples enhanced 
* signal (a,b), p^vr" signal (c,d) and p'V signal (e,f). 
The solid curve represents a projection of the maximum like- 
lihood fit result. The dashed curve represents the contribu- 
tion from continuum events, and the dotted line indicates the 
combined contributions from continuum events and B-related 
backgrounds. 



branching fractions to be 

► P+tt") = (10.9 ± 1.9 ± 1.9) X 10"^ 
>/?%+)= (9.5±1.1±0.8) X 10"^ 
^ /?%")= (1.4±0.6±0.3) X 10"^ 



B{B+ 
B{B+ 
BiB"^ 



where the first errors are statistical and the second sys- 
tematic. The systematic errors include the uncertainties 
in the efficiencies, which are dominated by the uncer- 
tainty in the tt^ reconstruction efficiency and in the case 
of p°7r+, by the uncertainty due to particle identification. 

Here we define the p^ir^ branching ratio by in- 

cluding those events that pass our selection and are fitted 
as signal but excluding those events that can be inter- 
j+TT^ with a p^, whose mass is closer 



preted as — 
to 0.77GeV/c^ than the mass of the reconstructed 
The signal significance for p°7r°, including statistical and 
systematic errors, is 2.1g, and we use a limit setting pro- 
cedure similar to Ref. [ll| to obtain a 90% Confidence- 
Level upper limit on its branching fraction. Fits on MC 
samples are used to find the signal hypothesis for which 
the ratio of the probablity that the fitted signal yield is 
less than that observed in data, and the probablity that 
the fitted yield is less than that in data under the null 
signal hypothesis, is 0.1. This signal hypothesis is shifted 
up by one sigma of the systematic error and the efficiency 
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is shifted down also by one sigma. This method gives an 
upper Umit of B{B° pV) < 2.9 x 10"^ 

Theoretical predictions of the ratio of branching frac- 
tions R = B{B° p^TT^)/B{B+ p°7r+), vary over 
a wide range. Tree level estimates suggest i? ~ 6 [l^ . 
while the inclusion of penguin contributions, off-shell B* 
excited states and scalar tt+tt" resonances leads to lower 
values, 7? ~ 2 — 3 0. Using the measured B^ — > p°7r+ 
branching fraction and the B'^ — )■ p^ir^ branching frac- 
tion from Ref. we find R = 2.38tg (stat)l°;^4(syst), 
which is in agreement with previous experimental re- 
sults 0. 

For the charged B decays we find the charge asymme- 
tries, A^^f = 0.24 ±0. 16 ± 0.06, yl^p"^ = -0.19±0.11± 
0.02, with contributions to the systematic errors listed in 
Table El 

In summary, we have presented measurements of 
branching fractions and CP-violating charge asymme- 
tries in 5+ p+tt" and _B+ p^ir^ decays, and a 
search for the decay 5° — > p'^ir^ . We observe the de- 
cay i?+ p+7r° with a statistical significance of 7.3f7. 
We also find a branching fraction for JB+ — > p'^ir^ that 
is consistent with previous measurements 0, and set an 
upper limit for B'^ p°7r°. We do not observe evidence 
for direct CP violation. 
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